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Abstract 

We investigate the parity-violating analyzing power in neutron capture on the proton 
at thermal energies in the framework of chiral effective field theory. By combining this 
analysis with a previous analysis of parity violation in proton-proton scattering, we are able 
to extract the size of the weak pion-nucleon coupling constant. The uncertainty is significant 
and dominated by the experimental error which is expected to be reduced soon. 



Although parity violation (PV) induced by the weak interaction is well understood at the 
level of elementary quarks, its manifestation at the hadronic and nuclear level is not that clear. 
This holds particularly true for the strangeness-conserving part of the weak interaction which 
induces PV in hadronic and nuclear systems. The Standard Model predicts PV forces between 
nucleons. However, their forms and strengths are masked by the nonperturbative nature of QCD 
at low energies. Combined with the difficulty of doing experiments with sufficient accuracy to 
extract parity-violating signals, hadronic PV is one of the least tested parts of the Standard 
Model. 

The understanding of low-energy strong interactions has increased tremendously by the use 
of effective held theories (EFTs). It has been realized that by writing down the most general 
interactions among the low-energy degrees of freedom that are consistent with the symmetries 
of QCD, one obtains an EFT, chiral perturbation theory (%PT), that is a low-energy equivalent 
of QCD. Each interaction in the chiral Lagrangian comes with a coupling strength, or low- 
energy constant (LEC), which needs to be extracted from data or computed in lattice QCD. In 
contrast to low-energy QCD itself, yPT allows one to calculate observables in a perturbative 
framework with expansion parameter pjA x , where p is the momentum scale of the process and 
A x ~ 1 GeV, the scale where the EFT breaks down. Although nuclear physics is intrinsically 
nonperturbative, the nucleon-nucleon (AW) potential can be calculated perturbatively within 
yPT. The resulting chiral potential is then iterated to all orders to calculate AW-scattering and 
bound state properties. This framework is usually called chiral nuclear EFT (for recent reviews, 
see Refs. M). 

The success of chiral EFT in parity-conserving (PC) nuclear physics has led to an analogous 
program in the PV sector |3 01013 Cl- One starts with the four-quark operators that are induced 
when the heavy weak gauge bosons are integrated out. The next step entails constructing a PV 
chiral Lagrangian which contains all interaction terms that transform under chiral symmetry in 
the same way as the underlying four-quark operators. From the resulting chiral Lagrangian one 
then calculates the PV AW potential and electromagnetic current. In the final step the obtained 
PV potential and current are applied, in combination with the PC chiral potential and current, 
in calculations of nuclear processes. The PV LECs appearing in the PV chiral Lagrangian can 
be fitted to some data and other PV processes can then be predicted. 

Although this sounds like a good strategy, in practice this procedure is complicated by the lack 
of data on PV processes. So far, hadronic PV has only been measured in a handful of experiments 
(see Refs. 00 for recent reviews). The longitudinal analyzing power (LAP), which would be 
zero in the limit of no PV, has been measured for proton-proton scattering at three different 
energies [HU HD EH, for proton-alpha scattering only at a single energy mm, and recently 
for the first time a preliminary result has been reported for radiative neutron capture on the 
proton np —> dy at thermal energies m- Nonzero parity-violating signals have also been found 
in more complex systems, as exemplified by the radiative decay of the 19 F nucleus [I6l [T? j and 
the anapole moment of the Cesium atom [18] . 

The first full chiral EFT analysis of PV nuclear forces has been done in Ref. [3] where it has 
been concluded that at leading order (LO) only a single interaction term appears: 

C ? = ^=W(7f x f) 3 W , (1) 

written in terms of the pion isospin-triplet 7r, the nucleon isospin-doublet N = ( p , n ) T , and the 
weak pion-nucleon coupling constant h n . The leading order PV potential arising from one-pion 
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exchange takes the form 
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in terms of the momentum transfer k = p — p' , where p and p' are the incoming and outgoing 
nucleon momenta in the center-of-mass frame, and <71,2 and Tip. the nucleon spin- and isospin- 
operators, respectively. F n = 92.4MeV denotes the pion decay constant, = 139.57MeV the 
charged pion mass, and = 1.29 the nucleon axial-vector coupling constant taking into account 
the Goldberger-Treiman discrepancy, in order to represent the strong 7 rA7V-coup 1 ing. 

Considering that there are no other terms at leading order, the one-pion exchange (OPE) 
potential can be expected to give the dominant contribution to PV in nuclear processes. Nev¬ 
ertheless, despite decades of experimental effort the existence of a long-range PV NN force has 
not been confirmed. This indicates that h n could be smaller than expected from naive dimen¬ 
sional analysis which predicts h 7r ~ GfF^A x ~ ICE 6 (consistent with the often-used estimate 
h n = 4.6 • ICE' of Ref. [19]), with Gf — 1-67 • ICE 5 GeV ~ 2 the Fermi coupling constant. In 
fact, the isovector nature of the weak pion-nucleon coupling already gives a natural suppression 
of sin 2 6 W ~ 1/4 [3:, 20], while a large- N c analysis indicates that h n is even further suppressed 
[201 EH [ 22 ]. A first lattice QCD calculation gave h n ~ 10 ' [23]. Finally, the absence of a PV 
signal in the 7 -ray emission from 18 F leads to the bound h n < 1.3 ■ ICE 7 )24[ [251126] . 

The evidence in favor of a small value of h n is not conclusive. Large-V c arguments can 
be misleading, especially for pionic interactions, while the lattice calculation did not include 
disconnected diagrams. The bound from 18 F depends on nuclear structure calculations of a 
relatively complicated nucleus and, despite being a careful work, might suffer from uncontrolled 
uncertainties. Finally, the Cesium anapole moment prefers a much larger value h n ~ ICE 6 
although the involved uncertainties are also larger [27, [28]. It seems that the only conclusive 
method of determining the size of h ^ is through a fit to experiments using simple few-body 
processes which are theoretically much better under control. Unfortunately, only a few PV 
signals have been measured so far in such few-body processes. In recent work we investigated 
the data on pp scattering in a chiral EFT framework El [29]. The main goal of this paper is to 
combine this analysis with the recent data on PV in radiative neutron capture on the proton 
ftp -A d'y and extract a value of h An analysis of PV in the inverse process yd —y np within 
pionless EFT has recently been performed in Ref. m- 

Our task gets complicated by two things. First of all, the OPE potential in Eq. ([2]) changes 
the total isospin and does not contribute to pp scattering. The three data points still carry 
information on the size of h^ because the analyzing power does depend on h n through the two- 
pion-exchange (TPE) diagrams [5[ [29 . 7]- The TPE diagrams appear at higher order in the 
chiral counting where additional contributions in the form of PV NN contact terms appear as 
well E1EHE2]- Secondly, although the PV OPE potential does contribute to np -A d'y capture, 
if the coupling constant h n is really as small as suggested, formally higher-order corrections 
can become relevant and need to be taken into account. Again such corrections appear as NN 
contact terms. We discuss these subleading terms in the PV potential and the current at a later 
stage. 

The other ingredients required for the calculation of PV observables are the PC NN potential 
and the PC and PV electromagnetic currents. As in Ref. [29] . we apply here the next-to-next-to- 
next-to-leading order (N 3 LO) chiral EFT potential obtained in Ref. [S3] and we refer the reader 
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to this paper for all further details. The N 3 LO potential exists for several values of the cut-off 
needed to regularize the scattering equation. Here, we regularize the PV potential in the same 
way as the PC potential via 

Vp V {p, p') e~ p6/A6 V PV (p, p')e~ p ' 6/A6 , (3) 


where three choices for A = {450, 550, 600} MeV are applied, see Ref. [iZSJUS]. TPE diagrams 
are regularized with a spectral cut-off A# = {500, 600, 700} MeV [33] . In recent work [33] 
an alternative regularization scheme (formulated in coordinate space) has been proposed which 
better preserves the long-range nature of pion-exchange terms in the potential. Considering the 
large experimental uncertainties in the field of nuclear parity violation, we do not expect drastic 
changes if the alternative regularization scheme is applied. Nevertheless, we will investigate this 
new scheme and its extension to N 4 LO [35] in future work. 

Within the chiral EFT power-counting rules the dominant PC current arises from the nucleon 
magnetic moments. At next-to-leading (NLO) order we encounter the one-body convection 
current)}] which arises from gauging the nucleon kinetic energy term, and the leading OPE 
two-body currents. The total PC current up to NLO is then given by 


Jpc = {~ [^s + /hvf] x 0 ) + (1 +Tj)(Pj +P/)}d (3) [P i -P- - q] 
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where fi s = 0.88 and / jl v = 4.72 are the isoscalar and isovector nucleon magnetic moments. 
The momenta of the incoming and outgoing nucleon interacting with the photon (of outgoing 
momentum q) are denoted by Pj and Pj, respectively. The momenta carried by the intermediate 
pions are k + q/2 = P\ — P{ and k — q/2 = P 2 7 — P 2 - In contrast, the leading PV current is solely 
due to OPE diagrams where one of the pion-nucleon vertices is from Eq. @ 
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These ingredients are sufficient to calculate the LO contribution to the longitudinal analyzing 
power in np —> dj. The details of the actual calculation will be presented in a longer paper m 
and therefore we focus here just on the results. 

1 Here the power-counting rules of Ref. m are followed where recoil and relativistic corrections are relegated 

to higher order by counting 1/mjv ~ fe/Aj, where k is the typical momentum scale of the process. The magnetic 

moment operator is not a recoil correction and only scales as 1/miv for conventional reasons. We thus treat 

Ps,v/mN ~ 1/A X which is also justified by the large value of p v = 4.72. 
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Table 1: Total cross section in mb for unpolarized np capture at 2.52 • 10 ~ 8 MeV lab energy. 
The first column is the calculated cross section using N 3 LO chiral potentials and the isovector 
nucleon magnetic moment p v . The second column also includes the leading PC OPE current. 
The experimental result is from Ref. |36j . Contributions from other currents at this order such 
as the isoscalar magnetic moment and convection current are negligible. 



isovector magnetic moment 

+ PC OPE currents 

Experimental result 

^tot 

305 ±4 

319 ±5 

334.2 ±0.5 


Table 2: Contributions to the LAP a 7 in np capture in units of h n . Part 1 is the contribution 
from one-body currents only, Part 2 from the isovector magnetic moment in combination with 
the PC OPE currents and the PV OPE potential, and Part 3 from the interference of the 
isovector magnetic moment and the PV OPE currents. 



Part 1 

Part 2 

Part 3 

Total 

Q"-y ! h'Jl 

-0.27 ±0.03 

-0.53 ±0.02 

0.72 ±0.03 

-0.11 ±0.05 


The longitudinal analyzing power (LAP) in np -± d r y capture is defined as 

. da + {6) - da_(0) 

AJ6) = -——- —— = a 7 cos 6 , 

7 da+(0) + da-(9) 7 


( 6 ) 


with da±(9) the differential cross section for incoming neutrons with positive/negative helicity 
and 9 the angle between the outgoing photon momentum q and the neutron spin. The experiment 
takes place at thermal energies where the total cross section for np capture is dominated by 
the nucleon isovector magnetic moment fi v , while the isoscalar magnetic moment p s and the 
convection current give negligible contributions. At NLO the PC OPE currents add to the cross 
section at the 5% level as can be seen in Table [I] For comparison, using the AV18 interaction we 
obtain 324 mbar in good agreement with Refs. [3HI 139] . The remaining discrepancy of roughly 
4% with respect to the experimental result should be removed by higher-order corrections, for 
example in the form of PC contact and TPE currents. In phenomenological models indeed 
the remaining discrepancy is explained by heavy-meson-exchange currents [38]. The theoretical 
uncertainties (l%- 2 %) quoted in the table are obtained from varying the cut-off parameters in 
the N 3 LO potential. 

Even in the presence of the PV potential, the numerator in Eq. (| 6 j) vanishes if we only 
include the leading magnetic moment currents. An interference with electric dipole currents, 
which appear at NLO in the form of the convection and OPE currents, is necessary to obtain a 
non-vanishing result. The dominant contributions to a 7 then consist of an interference between 
the isovector magnetic moment and: 


1. the one-body convection current in combination with the PV OPE potential, 

2. the two-body PC OPE currents in combination with the PV OPE potential, 

3. the two-body PV OPE currents. 
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All these contributions appear at the same order in the chiral counting and we present the re¬ 
sults in Table [2} The individual contributions are all of the same order, as expected from the 
power counting, and suffer only from minor uncertainties due to cut-off variations. However, the 
total result has a much larger relative uncertainty due to cancellations between the individual 
contributions. These cancellations were found also in Ref. uni where the AV18 potential has 
been applied in combination with the same currents. Our central value is also in good agreement 
with results based on various phenomenological strong potentials and the Siegert theorem for 
the electric dipole currents pumsg. These calculations do, however, not provide an uncer¬ 
tainty estimate. In Ref. |43| a smaller uncertainty was found when varying the cut-off (roughly 
±0.015 h n ), but the authors did not vary the strong potential simultaneously. In addition the 
Siegert theorem was applied for the electric dipole currents. The significant dependence of the 
total result on the cut-off parameters indicates that the extraction of h n from data on a 7 is less 
clean than might be expected. 

Having available the prediction for a 7 as a function of h^, we can now compare to data. For 
a long time only a bound on a 7 existed 

a 7 = (0.6 ± 2.1) • 10 -7 , a 7 = (-1.2 ± 1.9 ± 0.2) • 10 -7 , (7) 

from Refs. [41] and (55], respectively. Applying our most conservative estimate we obtain an 
upper bound \h n \ < 4.5 • 10 -6 on the weak pion-nucleon coupling. Recently a first preliminary 
result for a 7 was reported [T5] 

a 7 = (-7.14 ±4.4) • 10 ~ 8 . ( 8 ) 

This result is based on a subset of the full data taken and an improved result with an uncertainty 
at the 10~ 8 level is expected in the near future. 

To fit h n we combine the a 7 analysis with that of the LAP in pp scattering. As mentioned, 
the pp LAP does not depend on the OPE PV potential in Eq. ([2]) due to its isospin-changing 
nature. Nevertheless, the pp LAP still depends on h n due to the TPE potential which appears 
at NLO. At the same order the PV potential contains PV NN contact terms, but only one 
combination with LEC C contributes to pp scattering (C = —Co ± C\ ± C 2 — C 3 in the notation 
of Ref. [ 6 ]). The LECs h n and C were fitted to the pp data in Refs. [29]. In Ref. [ 6 ] the fit was 
slightly improved by including the dominant piece of the N 2 LO PV potential which does not 
depend on additional unknown LECs. 

We compare the obtained values of h n and C from the data on the np and pp LAPs in Fig. [l] 
using the intermediate cut-off values to regularize the potential. In the left panel, the ellipse 
denotes the contours of a total x 2 = 2.71 corresponding to a fit to pp data, whereas the red 
vertical lines denote a fit to a 7 at the one-sigma level. In the right panel, the ellipses denote 
contours of a total x 2 = 2,3,4 corresponding to a combined fit to the pp data and a 7 . The 
dashed lines in both panels are obtained if we use the expected future experimental uncertainty, 
±1 • 10 -8 , of the measurement of a 7 while using the same central value as in Eq. Q . The dashed 
contours are only there to illustrate what the accuracy could be with better data but should not 
be used to extract values of the LECs. 

From the right panel of Fig. [l] at the level of a total x 2 = 4 the LECs h n and C become 

K = (0.80 ± 0.70) • 10“ 6 , C = (- 6.0 ±3.0) • 1 (T 6 . (9) 

Since a x 2 analysis with so few data points can be misleading we collect in Table[3]the observables 
for the pp and np systems using three different fit values that all lie within the contours of 
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Figure 1: The allowed ranges for the LECs h^ and C (both in units of 10“ 6 ). Left panel: The 
blue ellipse is a fit to the pp data with a total % 2 = 2.71 and the vertical solid lines the fit of 
h n to a 7 at the one-sigma level. The vertical dashed lines correspond to the same fit, but now 
using the expected experimental uncertainty (±1 • 1CT 8 ) keeping the central value as in Eq. ([8]) . 
Right panel: The black (solid) ellipses are fits to the combined pp and np data with a total 
X 2 = 2,3,4. The red (dashed) ellipses are the same, but using the expected future experimental 
uncertainty. 


Fig. 0 From Table [ 3 ] we see that small (large) values of h n = 10 7 ( h ^ = 1.5 • 10 6 ) still give a 
reasonable fit to the pp data, but underpredict (overpredict) a 7 . However, considering the large 
experimental uncertainty of a 7 the fits cannot be excluded at a significant level. 

To study the dependence of the extraction of the LECs on the details of the strong NN poten¬ 
tial, we repeat the analysis for other cut-off values. For A = {450, 550, 600} MeV, respectively, 
the following best fit parameters emerge 

K = {0.48, 0.77, 1.1} • 10“ 6 , 

C = -{4.3, 6.4, 7.4} • 1(T 6 . (10) 

We see that the uncertainty in h ^ due to cut-off variations is roughly 40%. At present, the 
experimental uncertainty is still larger. However, once the precision of the a 7 measurement 
is improved by roughly a factor 2 the cut-off dependence will dominate the uncertainty. The 
observed cutoff dependence in the predicted value of a 7 is likely to be considerably reduced by the 
inclusion of higher-order corrections to the exchange currents. One frequently used approach 
along this line is to make use of the Siegert theorem. We emphasize, however, that such an 
approach yields only incomplete results for the exchange currents. We postpone a detailed 
investigation of the role played by higher-order contributions to a future study. 

In Fig. [ 2 ] we show contours for a total y 2 = 4 for three different cut-off values. To obtain an 
allowed range of the LECs, we extract the minimal and maximal values allowed by the three 
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Figure 2: The allowed ranges for the LECs h n and C (both in units of 10” 6 ). The ellipses 
correspond to fits to the combined pp and np data with a total x 2 = 4 for the three different 
cut-offs applied (from the smallest to the largest ellipse the cut-off changes from 450 to 600 
MeV). 


contours. With this conservative approach we obtain the following ranges 

K = { 1.1 ± 1 . 0 ) • 10” 6 , C = (-6.5 ±4.5) • 10” 6 . ( 11 ) 


The fits indicate that small values of h n ~ 10”' are barely consistent with the data, with values 
of h n ~ (5 — 10) ■ 10 ' being preferred. Such larger values disagree with the upper limit from 
1S F gamma-ray emission h n < 1.3 • 10”' and lattice and model calculations of h n — 10”'. An 
increase in the accuracy of the a 7 measurement is needed to make a firmer statement about this 
discrepancy. 

So far, our analysis included only the LO contribution to a 7 proportional to h^. If for 
whatever reason h n is small, formally subleading contributions might actually be dominant. 
The first corrections to a 7 appear two orders down in the chiral expansion and in principle 
consist of two contributions. The first arises from TPE diagrams in both the PV potential and 
currents. However, these contributions are proportional to h n as well and are then additionally 
suppressed by the assumed smallness of h n (in Ref. [43] TPE contributions were found at the 
10% level with respect to the OPE result based on the Siegert theorem and a phenomenological 
NN model). The other corrections appear in the form of PV NN interactions which contribute 
both to the PV potential and current. It turns out that the relevant potential and current 
depend on the same LEC C 4 , in the notation of Ref. [ 6 ], which is independent of the LEC C 
appearing in pp scattering. 


Vpy.NLO 

</pv,NLO 


—X t 2 ) 3 (<ti + <r 2 ) • k , 

~YJ2 ( fl • “ T 1 T 2 3 ) (<R + *2) . 


( 12 ) 

(13) 
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Table 3: Predictions for the pp LAP A z and a 7 (both in units of ICR') for three fits using the 
intermediate cut-off combination. The first fit corresponds to the best-fit value with h n = 0.77 
and C = —6.4. The second and third fits correspond to the values at the edge of the contours 
with h n = 0.1 and C = —3.0, and h n = 1.5 and C = —9.0 respectively (all in units of 10~ 6 ). The 
first three columns correspond to the pp LAP at three different energies and the fourth column 
to a-y. The experimental results are from Refs. mmmm- 



A z (13.6MeV) 

A z (45 MeV) 

A z (221 MeV) 

a 7 

Fit 1 

-0.90 

-1.56 

0.57 

- 0.74 

Fit 2 

-0.65 

-1.36 

0.50 

- 0.10 

Fit 3 

-0.89 

-1.19 

0.43 

- 1.44 

Exp. 

-0.93 ±0.21 

-1.50 ±0.22 

0.84 ±0.34 

-0.71 ±0.44 


Other contact contributions to the PV potential and current appearing at this order give rise to 
negligible contributions to a 7 . 

We obtain the total result for the asymmetry 


= (-0.11 ± 0.05)/i* + (0.055 ± 0.025)C 4 


(14) 


To estimate the size of the C 4 contributions to a 7 we can use resonance saturation. By compar¬ 
ison with the meson-exchange model of Ref. [T9], usually called the DDH model, the LEC C 4 
can be expressed ai 0 


F A 2 

C 4 = 

2 mjv 


ffcohi , gp(hy-h^) 


± 


mf 


mf 


(15) 


in terms of the masses m p — m u — 780 MeV and PC couplings g u = 8.4 and g p = 2.8. We have 


checked that both the potential and current in Eq. (12) m depend on this combination of DDH 
parameters by comparing to the currents derived in Ref. m- The sizes of the PV couplings 
h4,, hp, and h} p ' are unknown but can be estimated, albeit with significant uncertainty. Taking 
into account the whole reasonable range for these couplings as obtained in Ref. |19j . we find 
C 4 = (—0.8 ± 0.4) • 10 '. This range includes the more accurate prediction C 4 = — 1.2 • 10~' of 


Refs. [211[22l ]. To be conservative, we insert the DDH range into Eq. (14) to obtain the estimated 
uncertainty due to the short-range PV NN interaction 


= (-0.11 ± 0M)K - (0.5 ± 0.5) • 10" 


(16) 


Considering the current experimental uncertainty of ±4.4 • 10 -8 , the contact terms provide only 
a minor error. This would imply that the above analysis and extraction of h n is reliable. Small 
values of h n are thus disfavored, although formally not (yet) inconsistent. An improvement 
in the measurement of a 7 will provide a more definite answer regarding the size of h n . In 
Refs. [39i H2, 45j the dependence of the asymmetry on the short-range DDH parameters is found 


to be smaller than the central value of Eq. (14) by roughly a factor 4 to 5. A possible explanation 


2 In Ref. ( 6 | the resonance-saturation estimate of D 4 was found to also depend on h n due to TPE diagrams. 
However, in the calculation of a 7 we do not include TPE contributions explicitly so these terms should not be 
subtracted from the estimate. 

















Figure 3: The allowed ranges for the LECs C and C4 (both in units of 10 -6 ). The blue (solid) 
ellipse is a fit to the pp and np data with a total x 2 = 2.71 using h n = 0. The red (dashed) 
ellipse is similar but now using h n = 7.7 • 10“' which is the best-fit value for the intermediate 
cut-off. The dotted horizontal lines mark the predictions made by resonance saturation. 


might be the use of phenomenological strong potentials that typically have a stronger short- 
range repulsion than the chiral EFT potential, leading to a smaller dependence on short-range 
operators. A similar effect was found in the study of electric dipole moments |46| . 

The above reasoning is to some extent circular. We want to fit the LECs from few-body 
data only, but to estimate the effects of the formally subleading correction we require a model 
estimate of C 4 . This unfortunate situation is due to a lack of data which implies that we cannot 
fit all LECs at the same time in a consistent way. Some more insight can be obtained by using an 
alternative strategy. We force h n to be small and fit the LECs C and C4 to the pp and np data. 
We then obtain the fits in Fig. [3j The blue (solid) contour corresponds to a fit in the C — C 4 
plane with a total x 2 = 2.71 where we set h n = 0. In this case the fit prefers values for C4 which 
lie outside the range obtained from resonance saturation. The red (dashed) contour corresponds 
to a fit using h n = 7.7 • 10 -7 which corresponds to the best fit value for the intermediate cut-off 
combination. In this case the fit for C4 is centered around the zero and includes the resonance- 
saturation range marked by the dotted lines. Although strong conclusion cannot be drawn from 
this observation, it does indicate that small values of h n requires short-range contributions that 
are larger than expected. As always, more and/or more precise data are required to draw firmer 
conclusions. 

To summarize, in this paper we have extracted the values of two low-energy constants h n and 
C appearing in the parity-violating nucleon-nucleon potential and currents. To do so, we have 
used data on parity violation in proton-proton scattering and radiative neutron capture on a 
proton target. The extraction has been performed in the framework of chiral effective field theory 
which has been systematically applied to both the parity-conserving and parity-violating parts of 
the problem. We have estimated the uncertainties of the fits due to experimental uncertainties, 
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variation of cut-off parameters, and higher-order corrections and find the first of these to be 
dominant. Our extraction of the weak pion-nucleon coupling constant, h v = (1.1 ± 1.0) • 10~ 6 , is 
marginally consistent with bounds obtained in experiments on 18 F and a lattice QCD calculation. 
The expected increase in sensitivity of the a 7 measurement will significantly improve the fit and 
tell whether small values of h n are consistent with few-body experiments. 
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